str(fire_dataset)
## 'data.frame':    36011 obs. of  15 variables:
##  $ latitude  : num  -11.8 -11.8 -12.8 -14.4 -12.5 ...
##  $ longitude : num  142 142 133 143 131 ...
##  $ brightness: num  313 319 312 310 310 ...
##  $ scan      : num  1 1 3.1 1.1 4 1 1.2 1.2 1 3.4 ...
##  $ track     : num  1 1 1.7 1.1 1.9 1 1.1 1.1 1 1.7 ...
##  $ acq_date  : chr  "2019-08-01" "2019-08-01" "2019-08-01" "2019-08-01" ...
##  $ acq_time  : int  56 56 57 57 57 57 57 57 57 57 ...
##  $ satellite : chr  "Terra" "Terra" "Terra" "Terra" ...
##  $ instrument: chr  "MODIS" "MODIS" "MODIS" "MODIS" ...
##  $ confidence: int  48 71 42 33 36 68 24 79 72 40 ...
##  $ version   : num  6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 ...
##  $ bright_t31: num  297 297 299 296 299 ...
##  $ frp       : num  6.6 11.3 23.1 6.5 27.6 9.3 5.9 20.4 9.9 27.3 ...
##  $ daynight  : chr  "D" "D" "D" "D" ...
##  $ type      : int  0 0 0 0 0 0 0 0 0 0 ...

mean_brg median_brg min_brg max_brg
328.7507 326.4 300 504.4

Transform as date

str(fire)
## 'data.frame':    36011 obs. of  15 variables:
##  $ latitude  : num  -11.8 -11.8 -12.8 -14.4 -12.5 ...
##  $ longitude : num  142 142 133 143 131 ...
##  $ brightness: num  313 319 312 310 310 ...
##  $ scan      : num  1 1 3.1 1.1 4 1 1.2 1.2 1 3.4 ...
##  $ track     : num  1 1 1.7 1.1 1.9 1 1.1 1.1 1 1.7 ...
##  $ acq_date  : Date, format: "2019-08-01" "2019-08-01" ...
##  $ acq_time  : int  56 56 57 57 57 57 57 57 57 57 ...
##  $ satellite : Factor w/ 2 levels "Aqua","Terra": 2 2 2 2 2 2 2 2 2 2 ...
##  $ instrument: Factor w/ 1 level "MODIS": 1 1 1 1 1 1 1 1 1 1 ...
##  $ confidence: int  48 71 42 33 36 68 24 79 72 40 ...
##  $ version   : num  6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 ...
##  $ bright_t31: num  297 297 299 296 299 ...
##  $ frp       : num  6.6 11.3 23.1 6.5 27.6 9.3 5.9 20.4 9.9 27.3 ...
##  $ daynight  : Factor w/ 2 levels "D","N": 1 1 1 1 1 1 1 1 1 1 ...
##  $ type      : int  0 0 0 0 0 0 0 0 0 0 ...

no missing values

str(fire)
## 'data.frame':    36011 obs. of  15 variables:
##  $ latitude  : num  -11.8 -11.8 -12.8 -14.4 -12.5 ...
##  $ longitude : num  142 142 133 143 131 ...
##  $ brightness: num  313 319 312 310 310 ...
##  $ scan      : num  1 1 3.1 1.1 4 1 1.2 1.2 1 3.4 ...
##  $ track     : num  1 1 1.7 1.1 1.9 1 1.1 1.1 1 1.7 ...
##  $ acq_date  : Date, format: "2019-08-01" "2019-08-01" ...
##  $ acq_time  : int  56 56 57 57 57 57 57 57 57 57 ...
##  $ satellite : Factor w/ 2 levels "Aqua","Terra": 2 2 2 2 2 2 2 2 2 2 ...
##  $ instrument: Factor w/ 1 level "MODIS": 1 1 1 1 1 1 1 1 1 1 ...
##  $ confidence: int  48 71 42 33 36 68 24 79 72 40 ...
##  $ version   : num  6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 ...
##  $ bright_t31: num  297 297 299 296 299 ...
##  $ frp       : num  6.6 11.3 23.1 6.5 27.6 9.3 5.9 20.4 9.9 27.3 ...
##  $ daynight  : Factor w/ 2 levels "D","N": 1 1 1 1 1 1 1 1 1 1 ...
##  $ type      : int  0 0 0 0 0 0 0 0 0 0 ...
ggplot(fire, aes(x = acq_date, y = brightness))+
  geom_point(aes(color = daynight, shape = daynight)) +
  scale_shape_manual(values=seq(0,1))

str(fire)
## 'data.frame':    36011 obs. of  15 variables:
##  $ latitude  : num  -11.8 -11.8 -12.8 -14.4 -12.5 ...
##  $ longitude : num  142 142 133 143 131 ...
##  $ brightness: num  313 319 312 310 310 ...
##  $ scan      : num  1 1 3.1 1.1 4 1 1.2 1.2 1 3.4 ...
##  $ track     : num  1 1 1.7 1.1 1.9 1 1.1 1.1 1 1.7 ...
##  $ acq_date  : Date, format: "2019-08-01" "2019-08-01" ...
##  $ acq_time  : int  56 56 57 57 57 57 57 57 57 57 ...
##  $ satellite : Factor w/ 2 levels "Aqua","Terra": 2 2 2 2 2 2 2 2 2 2 ...
##  $ instrument: Factor w/ 1 level "MODIS": 1 1 1 1 1 1 1 1 1 1 ...
##  $ confidence: int  48 71 42 33 36 68 24 79 72 40 ...
##  $ version   : num  6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 6.3 ...
##  $ bright_t31: num  297 297 299 296 299 ...
##  $ frp       : num  6.6 11.3 23.1 6.5 27.6 9.3 5.9 20.4 9.9 27.3 ...
##  $ daynight  : Factor w/ 2 levels "D","N": 1 1 1 1 1 1 1 1 1 1 ...
##  $ type      : int  0 0 0 0 0 0 0 0 0 0 ...
library(gganimate)
p_1 <- ggplot(
  fire, 
  aes(x = scan, y = brightness , colour = satellite)
  ) +
  geom_point(show.legend = FALSE, alpha = 0.7) +
  scale_color_viridis_d() +
  scale_size(range = c(0, 1)) +
  scale_x_log10() +
  labs(x = "Day or night ", y = "brightness")

p_1  + transition_time(acq_date) + ggtitle('Eathquake magnitude type defined per Magnitude and per year ',
          subtitle = 'acq_date: {frame_time}')